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We analyse the present experimental evidence for a complex CKM ma¬ 
trix, even allowing for New Physics contributions to ex, o-j/^Ks^ 

AMb^, and the AI = 1/2 piece oi B ^ pp and B pir. We emphasize the 
crucial role played by the angle 7 in both providing irrefutable evidence for 
a 3x3 complex CKM matrix and placing constraints on the size of NP con¬ 
tributions. It is shown that even if one allows for New Physics a real CKM 
matrix is excluded at a 99.92% C.L., and the probability for the phase 7 
to be in the interval [—170°;—10°] U [10°; 170°] is 99.7%. Large value of 
the phase y, e.g. of order A, is only possible in models where the unitarity 
of the 3x3 Cabibbo-Kobayashi-Maskawa matrix is violated through the 
introduction of extra (5 = 2/3 quarks. We study the allowed range for 
X and the effect of a large y on various low-energy observables, such as 
CP asymmetries in B meson decays. We also discuss the correlated effects 
which would be observable at high energy colliders, like decays t cZ, 
etc.. 


1. INTRODUCTION 

As is well-known Unitarity Triangle fits indicate the prominent role 
of the Cabbibo-Kobayashi-Maskawa (CKM) mechanismpp in CP violation 
and Flavour Physics. The huge amount and the variety of CP collected 
data allows for a systematic search of Physics Beyond the Standard Model 

(SM)[2IS1- 

A quite general and natural framework to go beyond the SM is: 

1. Allow for New Physics in every place except in weak tree-level domi¬ 
nated processes. 

2. Assume 3x3 unitarity. 

This framework is general enough to include practically all models of 
New Physics (NP) except those that explicitly violate 3x3 unitarity. Going 
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beyond 3x3 unitarity makes the analysis almost impossible in a model 
independent way|l]. We will proceed in the first part with assumptions 1 
and 2 . 

To go on with this analysis it is very important to understand the key 
role played by a measurement of the unitarity triangle phase 7 . It is very 
convenient to write the CKM matrix in the following way|n][S]: 

/ \Vud\ \Vus\e^^' \Vub\e-^^ ... \ 

-\Vcd\ IW.I iWfel ... 

\Vtd\e-^f^ \Vtb\ ... 

V . 

Where the rephasing invariant CP violating phases arejZj 

P = avg{-V,,V:,V;,Vt,) , j = ^Tg{-V^dV:,V:,Wt) .21 

X = s^rg{-VtsV*V:,V,b) , x' = arg(-WdK:K‘dK.) 

It has been shown from 3x3 unitarity -A ~ 0 . 2 - that x' ~ and X ~ 
even outside the SM|H1. So x' is too small to consider its measurement. 
(3 and x, accompanying \Vtd\ and |I4s| respectively, will enter only in loop 
processes: virtual transitions q ^ t. Therefore processes that in the SM 
measure (5 and x could be contaminated by NP. Processes that measure 
f3 in the SM framework will be no longer evidence of CP violation in the 
CKM matrix in the presence of NP. The unique measurable phase that can 
be extracted from tree-level processes and therefore not contaminated by 
NP is 7 , because it can appear in a tree-level transition 6 ^ n j^J. 

We can conclude that just 7 and the moduli of the first two rows are the 
unique parameters whose extraction from experimental data is not contami¬ 
nated by NP. Or put in another way, these parameters -extracted from weak 
tree-level decays- are valid in all the models included in our assumptions 1 
and 2 . 

In section II we will parametrize and clarify the NP physics contributions 
to the contaminated flavour and CP observables. Having seven experimental 
data of the CKM matrix in these general class of models, in section III 
we will address the question of the dominance of the CKM mechanism 
as the origin of CP violation in these models. At the same time, with 
the contaminated observables, we will set bounds on the NP contributions. 
Finally in section IV we will study potential large deviations of x from its 
SM value in models that violates 3x3 unitarity. 

2. NP CONTRIBUTIONS IN a, (3 AND 7 MEASUREMENTS 

As we have explained the measurements of (3 will be contaminated by 
NP. Also it is known that (3 measurements come always from the f3 con- 
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tribution to — B^d mixing. The simple reason why (5 contributions to 
B decay amplitudes cannot be measured - in a clean theoretical way- is 
because it always enter in the SM with a second weak phase. Therefore it 
is enough to parametrize the NP contributions to the B^ — B^ mixing. We 
follow reference jH] to write the off diagonal mixing matrix as 

(3) 

the SM corresponds to = 1 and = 0. In this way we get 


q 

P 


Bd 



SM 

Bd 


= g-2*(/3-0d) 


( 4 ) 


and therefore, by defining 




13 = fi-(t^d 

( 5 ) 

we have for the experimental 

observables 


AMs, 

= 

( 6 ) 

Sj/pKs 

= sin 2 (/3 — (^ci) = sin 2/3 

( 7 ) 


Clearly an independent knowledge of the full CKM matrix together with 
the B'^ —> J/'il’Ks asymmetry Sjj^Xs the mass difference will 

give us the opportunity to test for the NP parameters Vd and cpd- 

As we will see in the next section, the actual knowledge of the moduli 
of the CKM matrix in the first two rows is not enough to know the entire 
CKM matrix, so the measurement of gamma is of paramount importance in 
order to complete the knowledge of the CKM matrix and to proceed with 
the NP analysis. Gamma can be obtained from the phase of the rephasing 
invariant quartet VusVcbV*fjV*g (H]. A way of measuring it from pure tree- 
level decays is trough the interference of the two pure charged currents decay 
paths b —> sue and b scu ^njjTJ. Babar^2] and Belle^Sl have presented 
results using the Dalitz plot analysis in B^ —> DK^ with the subsequent 
decay D^,D —> KsTT~^'7r~. If any potential NP in — D mixing is 
neglected|14|. as seems reasonable, these analysis provide a measurement of 
gamma free from NP contributions. 

Other relevant ways of measuring 7 are the methods to measure 
a = arg(-K,,K;,K;,Kfe) M- By now it is well-known that using eq. m 
one has a = 7 r — /3 — yby definition|lfij|17j. So a measurement of a is noth¬ 
ing else than a measurement of /3 -|- 7 - The main channels are B tttt, pvr 
and pp. In this case, the presence of penguin pollution could become NP 
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pollution, so one has to be much more careful in these processes. The rele¬ 
vant observables are 


^fq\ ^ /) 


( 8 ) 


To understand the effects of NP let us first neglect the SM penguin pollution 
and treat the tt-tt channel as if it were a pure tree level and therefore without 
NP pollution in the decay amplitudes. In these case we have in the 7r'’'7r“ 
channel 


A+_ = e“*2/3e-*27 = e*2a 


(9) 


where a = tt — (3 — ^ = a + (fid- The CP asymmetry measures a. Once we 
include penguin pollution, the Gronau and London isospin analysis can be 
partially summarized with the following formula^ [TK] 


A+o = - 


A{B- 


vr vr^) 


1 R ± iy — K 




( 10 ) 


where 


<-i Kop+^iALi'-iiooi" 

^ =-cl- 

y = i^+oiiAi+-i (12) 

that tell us that A+o can be extracted from the experimental data including 

branching ratios and A^_ But A^q in the SM and therefore in our NP 

scenario is a pure tree-level amplitude with weak phase gamma - is propor¬ 
tional to the AI = 3/2 piece-, so we will have that the observable A+o will 
be _ 

A+o = e-*2/3e-*27 = e-2ia 

We conclude that the usual way of measuring a in vrvr decays provides 
us with a measurement of a = tt — /9 — 7 even in the presence of NP 
in the AI = 1/2 Diece|19j. Obviously the knowledge of (5 from J/tpKs 
converts these a methods in another way of extracting 7 from tree-level 
pieces. Similar results are obtained for pvr and pp. 

3. IS THE CKM MATRIX COMPLEX IN THE PRESENCE 

OF NP? 

Within this class of models, in order to investigate whether the present 
experimental data already implies that CKM is complex, one has to check 
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whether any of the unitarity triangles is constrained by data to be non- 
“flat”, i.e. to have a non-vanishing area. If any one of the triangles does not 
collapse to a line, no other triangle will collapse, due to the remarkable prop¬ 
erty that all the unitarity triangles have the same area. This property simply 
follows from unitarity of the 3x3 CKM matrix. The universal area of the 
unitarity triangles gives a measurement of the strength of CP violation medi¬ 
ated by a VT-interaction and can be obtained from four independent moduli 
of VcKM- The fact that one can infer about CP violation from the knowl¬ 
edge of CP-conserving quantities should not come as a surprise pH] . It just 
reflects the fact that the strength of CP violation is given by the imaginary 
part of a rephasing invariant quartet pi] , J = ±Im with 

{i ^ j,a ^ (3), which in turn can be expressed in terms of moduli, thanks 
to 3x3 unitarity. Restricting ourselves to the first two rows of Vckm, to 
avoid any contamination from NP, a possible choice of independent moduli 
would be I Pus I, iPcfel) |Pufe| and \Vcd\- One can then use unitarity of the first 
two rows to evaluate J, which is given, in terms of the input moduli[2j, by 


4J2 = 4 (l - iPufeP - iPusP) iPubP iPcdl" iKfeP - 

- {\Vusf - iPcdl' + iKdP \Vubf - |Pc6p \Vubf - |Pc6p |P.sP)a4) 

Note that Eq. m is exact, but the actual extraction of J from the chosen 
input moduli, although possible in principle, it is not feasible in “practice”. 

To illustrate this point, let us consider the present experimental values of 
I Pus I ) iPcftI ) and |Pcd|) assuming Gaussian probability density distribu¬ 
tions around the central values. We plot in Fig. o the probability density 
distribution of J^, generated using a toy Monte Carlo calculation. Only 
31.1% of the generated points satisfy the trivial normalization constraints 
and, among those, only 7.9% satisfy the condition that the unitarity trian¬ 
gles close (> O). 

As we have mentioned, including the experimental data Sjj^Ks 

still does not give evidence of a complex Vckm ■ We cannot con¬ 
clude if Vckm is the dominant contribution to CP violation, although we 
can set bounds on the NP parameters and (l)d- In Fig.( 2(a) I we plot 68% 


(black), 90% (dark grey) and 95% (grey) probability regions of the probabil¬ 
ity density function (PDF) of the apex —VudV^^/VcdV*^ of the db unitarity 
triangle. In Fig.( |2(b)] ) we represent joint PDF regions in the plane (r^, 2(pd)- 
Because 7 gives the apex of the triangle, it is clear from Fig.( |2(a)] ) that there 
is essentially no restriction on 7 . On the contrary, because the moduli of 
the first two rows put an upper bound on |/3| and u pper and lower bounds 
on Rt = iPdPjftl / iPcdPcftI) we can see in Fig.( |2(b)1 ) significant constraints 
on 2(j)d and r^. 
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Fig. 1. J2 distribution from |Ks| = 0.2200 ± 0.0026, |Kb| = (3.67 ± 0.47) x lO'^, 
\Vcb\ = (4.13 ± 0.15) X 10-2 and |Kd| = 0.224 ± 0.012. 


E 


T3 

CN 


(a) Apex of the unitarity triangle (b) {^r^,2(j)d) joint distribution. 

bd. 

Fig. 2. Probability distributions, no restriction on 7 . 




Including the measurement of gamma in Fig.( |2(a)] ) will fix the unitarity 
triangle and the Vckm matrix. The enormous effort developed at the B- 
factories Belle and BaBar has resulted in the first measurements of 7 in tree- 
level decays —> DK^, D*K^ -^^5 where the two paths 
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to or interfere in the common decay channel D^,D^ —> 

From a Dalitz-plot analysis , Belle has presented 7 = 68°± zb 13°±11° 
and BaBar ^21 7 = 70°±26°±10°±10°, together with the solutions obtained 
by changing 7 ^ 7 ± tt. 

We average conservatively both measurements to the value 7 = 69° ±21° 
(—111° ± 21°), which we take as a quantitative measurement of a complex 
CKM matrix independent of the presence of NP at the one-loop weak level. 

BaBar has also presented a time-dependent analysis of the p~^p~ chan¬ 
nel m, that once supplemented with the p^p^ and p^p^ branching ratios 
[23 mi and the measurement of the final polarization, can be translated 
[23 into the measured value a = 96° ± 10° ± 5° ± 11° where the last er¬ 
ror comes from the usual SU{2) isospin bounds^ la^ff— a\ < 11°. Be¬ 
cause the measurement is sensitive to sin(2ae//)) cf = cie// 11° presents 
a fourfold ambiguity (a, a ± vr, ^ — a, —a — f) • In the p-ir channel the pen¬ 
tagon isospin analysis - from quasi-two-body decays - needs more statistics 
and/or additional assumptions. A time-dependent Dalitz-plot analysis in 
the channel B 7r'*'7r“7r^ has been presented by BaBar m, with the re¬ 
sult a = 113° ± ± 6°.Since this analysis is sensitive to both sin {2aeff) 

and cos (2ae//), the resulting ambiguity is just a twofold one (a, a ± vr). It 
is remarkable that these two solutions are in good agreement with two of 
the solutions coming from the pp channel. This important property will be 
used to eliminate two of the four solutions coming from the pp channel. 

The situation in the Trvr channel does not yet allow a full isospin anal¬ 
ysis and the isospin bounds are quite poor. Furthermore BaBar and Belle 
measurements are still in some conflict, so that we will not use these results. 

As before, we average the data from pp and pvr but only keep the two 
solutions consistent with the pn channel data. Our averaged values are 
a = 100° ± 16°, (-80° ± 16°). 

To analyze the implications for the dominance of the CKM mechanism 
for CP violation and the presence of NP, we add both measurements 7 and 
a to the previous analysis presented in Figs.( |2(a')] ) and ( 2(b)| ). 

In Fig.( |3(a)'] ) we represent the analogue of Fig.( |2(a) ). We conclude [2j 
that a real CKM matrix is excluded at a 99.92% C.L. and the probability 
of 7 E [10° ; 170° ] U [-170°; -10°] is 99.7%. 

In Fig( |3(b)| ) we can see three solutions in the (r^,2i/)d) plane with 
2((>d ~ 0°, —75° and —150°. In these plots one has in general four solutions 
corresponding to the two values of 7 (a) and to the two signs of cos (2/3) 
The last solution corresponds to cos (2/3) < 0. It is the inclusion of both 7 
and a constraints that almost eliminates the cos (2/3) < 0 solutions. The 

^ In our notation f is the usual a^f / but where we have introduced (3 instead of /3 

as the phase in the B^-B^ mixing. 













E 



(a) Apex of the unitarity triangle (b) (r^,2(/)d) joint distribution. 

bd. 

Fig. 3. Probability distributions including 7 and a experimental data. 


first solution is obviously the SM one, and the semileptonic asymmetry Asl 
- not included here - is starting to play an important role [77] in the exclusion 
of the solution 2(l)d 75°. 

It is important to stress that having an irrefutable piece of evidence for a 
complex CKM matrix, in a framework where the presence of NP is allowed, 
has profound implications for models of CP violation. In the particular 
case of models with spontaneous CP violation, a complex CKM matrix 
favours the class of models where, although Yukawa couplings are real, the 
vacuum phase responsible for spontaneous CP violation also generates CP 
violation in charged-current weak interactions. Conversely, the evidence for 
a complex CKM matrix, even allowing for the presence of NP, excludes the 
class of models with spontaneous CP violation and a real CKM matrix at 
a 99.92% C.L.. 


4. THE SIZE OF x = arg {-VtsV*^y^Vcb) AND DEVIATION 
FROM 3x3 UNITARITY 

Within the SM and any extension where Vsxs is unitary, like supersym¬ 
metric or multi Higgs doublet models, we have the relation [7] 


\Vuh\\Vus 

iKcfellK,,! 


sin(7 -b x' - X) 


smx = 


(15) 
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which shows that |x| ^ in any model where 3x3 CKM unitarity holds. 
In particular, within the SM one obtains at 90% CL 


0.015 < X < 0.022 (SM) (16) 

The only models in which x can be significantly larger than are those 
in which Vsxs is not unitary, what can only be achieved by enlarging the 
quark sector. The most simple way of doing this is with the introduction of 
new quark singlets ISHl^ .Quark singlets often arise in grand unified theories 
and models with extra dimensions at the electroweak scale |29] . They have 
both their left- and right-handed components transforming as singlets under 
SU{2)l, thus their addition to the SM particle content does not spoil the 
cancellation of triangle anomalies. In these models, the charged and neutral 
current terms of the Lagrangian in the mass eigenstate basis are 

Cw = -^ULl^VdLW+ + h.c. (17) 

Cz = {ulI^Xul - - 2slrJ^^) (18) 

Zcw 

where u = {u,c,t,T,...) and d = {d, s,b, B,...), V denotes the extended 
CKM matrix and X = U = V^V are hermitian matrices. X and 

U are not necessarily diagonal and thus flavour-changing neutral couplings 
(FCNC) exist at the tree level, although they are naturally suppressed by 
the ratio of the standard quark over the heavy singlet masses. Moreover, 
the diagonal Zqq couplings, which are given by the diagonal entries of X 
and U plus a charge-dependent term, are also modified. Within the SM 
Xuu — Xf^c — X^i — 1 , Xqq! — 0 for q ^ q ■, Udd — Ugs — Ufii) — 1 and 
Uqq' = 0 for g 7^ g'. The addition of up-type Q = 2/3 singlets modifies the 
first two of these equalities, while the addition of down-type Q = —1/3 ones 
modifies the last two. 

In models with a down quark singlet, from orthogonality of the second 
and third columns of V, one obtains the generalization of Eq. (USD 


sinx = 


\Vub\\Vus 

\Vcb\\Vcs\ 


sin(7 + x' -X) 


Im {Ubse-^^) 


(19) 


From the present bound on 6 —> si~^t , one obtainsjnn][S2 that at most 
\Ubs\ 10“^ ~ thus implying that in this class of models x cannot be 
significantly larger than A^. 

^ The addition of a sequential fourth generation is another possibility, but it is dis¬ 
favoured by two facts: (i) the experimental value of the oblique correction parameters 
only leave a small range for the masses of the new quarks; (ii) anomaly cancellation 
requires the introduction of a new lepton doublet, in which the new neutrino should 
be very heavy, in contrast with the small masses of the presently known neutrinos. 
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In model with an up quark singlet, from orthogonality of the second and 
third rows of V, one getsjB] 


sinx = 


Im Xct 
\Vcs\\Vts 


+ 0(A2) 


( 20 ) 


In contrast with models containing down-type singlets, where the size of 
all FCNC is very restricted by experiment, present limits on are rather 
weak. The most stringent one, \Xct\ < 0.41 with a 95% CL, is derived 
from the non-observation of single top production at LEP, in the process 
e+e“ —> tc and its charge conjugate. This bound does not presently provide 
an additional restriction on the size of x- In models with extra up singlets 
I Xct I can be of order A^, yielding x ~ A. From Eq. one derives some 
important phenomenological consequences. First, we observe that a sizeable 
X is associated to a FCNC X^ ~ 10“^, which leads to FCNC decays t —> 
cZ at rates observable at LHC. In addition, the modulus of X^ obeys the 
equality m 

|Xct|2 = (l-Xcc)(l-Xtt) (21) 

This relation shows that conditions for achieving Xct ~ 10“^ are to have a 
small deviation O(A^) of Xcc from unity (which is allowed by the measure¬ 
ment of Rc and pIT| i and a deviation of Xu ~ from unity of order 
A^. This deviation is only possible if the mass of new top quark T is bellow 
1 TeV, again testable at LHC. Finally the — D mass difference sets 
bounds on Xuc that for large x - (1 ~ ^cc) ~ 0('^^) - translate into bounds 
on (1 — Xuu) in such a way that deviations of unitarity in the hrst row are 
not observable. By the same token, large values of x will be correlated with 
an important contribution to — D mixing |S|. 


The shaded area in Fig. represents the allowed interval of x for a given 
Xu- Note that lower values of Xu are allowed for lighter T quarks [SJ. The 
Fig.® has been obtained incorporating all the relevant constraints: the 
correction to the oblique parameter AT, Rc and in the charm sector, 

— 7^ mixing and several bounds from rare K and B decays. 

Because x is the phase of Vts, important effects will appear in 6 —> s 
transition in the piece nevertheless the presence of the new T quark 

will introduce a contribution proportional to VtsV^^ and dependent on the 
mass mj' of this quark. For the CP asymmetry of —> cftKs we get|H] 

S^Ks = sm(2;3 + 2x) (22) 

where 

.. 1 ( 1 + /(^T, mt)VT,V^jVM \ 

2 ^\l +f{mT,mt)V*^VTs/V*^Vts) 


(23) 
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Fig. 4. Allowed interval of x (shaded area) as function of Xu 


gets contributions from x and the new T quark. /{rriT, rrit) is in general a 
complex function fixed by Wilson coefficients and matrix elements [S]. By 
scanning all the allowed range of parameters we get S(j,Kg G [0.57,0.93]. 

Less dependent on hadronic matrix elements is the contribution to — 
Bg mixing: 

i,j=t,T' 

(24) 


2iXeff 


e-2ix 


^ S{mt,mT)VTs^Tb 

S{mt,mt)VilVtb) 


S{mT,mT) / S{mt,mT) \‘^ f ^Ts^Tb \^\ 
S{mt,mt) \S{mt,mt)) \V^gVtb ) j 


(25) 


with K a constant factor and S the usual Inami-Lim box function |33j. In 
any channel without a weak phase in the decay amplitude, for example in 
the Bg —> Df Dj and ^|:([> channels, the time dependent CP asymmetry is 


~ ^Ae// (26) 

again Xeff is equal to x plus a T dependent contribution. The range of 
both contributions goes in opposite direction as ttit changes. 
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Fig. 5. Range of variation of the asymmetry in terms of mr 


In Fig. ©, we present the range of variation of Sj^+j^- in terms of m^. 
A potential spectacular departure from the SM could be seen at LCH. 

5. CONCLUSIONS 

The first measurements of 7 clearly points towards the CKM mechanism 
as the dominant source of CP violation in the quark sector, even in the 
presence of NP in all loops. A real CKM matrix is excluded at the 99.92% 
C.L. Including a data, there are still alive three solutions: the most robust is 
the SM, one with cos (2/3) < 0 is almost excluded and it remains another NP 
solution where the semileptonic asymmetry can do a relevant job. Future 
improvements of the data will be crucial to left just very small deviations 
of the SM. 

A large deviation of the SM value of the phase y is only possible in 
models that violates 3x3 unitarity. More precisely, in models with an up 
singlet quark, x can be order A. Moderates effects appear in —> cjjKs-, 
more spectacular effects can be present in CP asymmetries in the system. 
Correlated effects of this scenario would be: rare top decays t —> cZ at a 
rate observable at LHC, production of a heavy top T at LHC and important 
deviations from 1 of Xu ~ measurable at ILC. 
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